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Abstract. Recent observations of r-nuclei, and in partic- 
ular of Th, in ultra-metal poor stars revived the old idea 
that the Th cosmochronometry could provide an age es- 
timate of the oldest stars in the Galaxy, and therefore a 
lower limit to the age of the Galaxy. Unfortunately, some 
nuclear, astrophysics and observational uncertainties still 
affect the theoretical r-process models required to predict 
the original production of Th. The impact of these un- 
certainties on the prediction of the age of the Galaxy is 
analyzed and discussed. 
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1. Introduction 

One of the methods called for estimating the age of the 
Galaxy is based on the analysis of observed abundance 
of some long-lived radionuclides. The most studied cos- 
mochronometries involve 187 Re or the actinides 232 Th, 
235 U and 238 U. Though some promising results have re- 
cently been achieved in decreasing the uncertainties affect- 
ing the 187 Re cosmochronometry, the predictions based 
on the trans-actinides can still be regarded as relatively 
poor (Arnould & Takahashi 1999). Nevertheless, the re- 
cent observation of r-process elements, including Th, in 
ultra-metal-poor halo stars, such as CS 22892-052 or HD 
115444 (Sneden et al. 1996, 1998) has brought some re- 
newed excitement in the estimate of the age of the Galaxy 
on grounds of the Th cosmochronometry (Cowan et al. 
1997; Pfeiffer et al. 1998). With a metallicity as low as 
[Fe/H] = -3 and a composition enriched in some pure r- 
elements, these stars provide strong evidence that the pro- 
duction of heavy elements by the r-process already took 
place early in the history of the Galaxy. Moreover, the 
abundance pattern of the 15 r-elements heavier than Ba 
at the surface of CS 22892-052 (or the 9 elements in HD 
115444) shows a striking similarity with the solar system 
r-abundance distribution, leading to the tempting (though 
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hazardous) conclusion that the r-process mechanism is 
"unique" , i.e any astrophysical event producing r-elements 
gives rise to a solar-like abundance distribution. This con- 
clusion has been critically analyzed by Goriely & Arnould 
(1997) who showed that this assumption may be valid in- 
deed, but is by far not the only possible one, as the obser- 
vations in the limited 56 < Z < 76 range are equally com- 
patible with an abundance distribution that does not fit 
the solar one outside the observed domain. This ambigu- 
ity is assigned to the fact that the observed CS 22892-052 
pattern of abundances reflects primarily nuclear physics 
properties, and not one or another specificity of a blend 
of r-process events. This universality assumption is a fun- 
damental prerequisite to build a Th cosmochronometry 
upon the abundance analysis of metal-poor stars at the 
present time. In principle, it could be possible to derive 
the abundance of Th ingested in these metal-poor stars 
from theoretical extrapolations based on direct fits to the 
observed abundances. However, in practice, this exercise 
is affected by uncertainties even greater than when bas- 
ing the fits on the solar abundances, because of the re- 
stricted number of elements observed, the impossibility to 
distinguish isotopic ratios and the smaller precision in the 
abundance determination compared with the data avail- 
able in the solar system. In particular, Goriely & Arnould 
(1997) showed that the r-elements distribution at the sur- 
face of CS 22892-052 could be reproduced satisfactorily 
by a random superposition of canonical r-process events. 
In this case, the theoretical extrapolation to the actinide 
region based on parametric r-process models is simply 
meaningless. Nevertheless, future accurate observations of 
r-elements in ultra-metal-poor star could change this sit- 
uation. 

For this reason, we will consider in the present paper 
the universality assumption to be valid in order to analyze 
if, despite this difficulty, the recent accurate observation 
of Th at the surface of ultra-metal-poor stars can indeed 
provide a reliable estimate of the stellar age by comparing 
it with the universal r-abundance of Th. Such a procedure 
requires the estimate of the Th by r-process models, which 
are known to suffer from very many astrophysics and nu- 
clear physics problems, in spite of much recent theoretical 
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and experimental effort. In this respect, the Th problem is 
particularly acute, since with U, Th is the only naturally- 
occuring nuclide beyond 209 Bi, so that the estimate of Th 
production relies on extrapolation procedures based on fits 
to the solar (or stellar) r-abundance distribution. In Sect. 

2, a brief description of the adopted r-process models is 
given in relation to the Th cosmochronometry. In Sect. 

3, the various uncertainties affecting the Th production 
are studied and their impact on the estimate of the stel- 
lar age is analyzed. In Sect. 4, it is shown by comparing 
the solar fits to the stellar r-element distribution observed 
that future observations of Pb, Bi or U could put the Th 
cosmochronometry on safer grounds. 

2. Th cosmochronometry and the r-process 

Assuming that the whole r-abundance distribution ob- 
served in CS 22892-052 and HD115444 is essentially solar, 
it is straighforward to relate the star age T* to the Th 
abundances, 

(iL - (S). <» 

where r(Th) = 20.27 Gyr is the characteristic a-decay 
timescale of Th and the subscripts obs and r refer to 
the observed and universal r-process abundance ratios, 
respectively. As classically done, the Th abundance is 
here expressed relative to the spectroscopically relevant 
Eu r-dominant element. The recent accurate observa- 
tion of Th at the surface of CS 22892-052 amounts to 
log(Th/Eu) o6s = -0.70 ±0.08 (Sneden et al. 1996; Cowan 
et al. 1997). Th has also been observed at the surface of 
HD 115444, but its precise abundance remains to be de- 
termined (Pfeiffer et al. 1998). Assuming that a solar-like 
mix of the r-elements ingested in these halo stars origi- 
nates from a small number of nucleosynthetic events that 
took place just before the formation of the stars, the age 
of the star can be estimated from Eq. ([!]) without call- 
ing for a complex model of the chemical evolution of the 
Galaxy. The only difficulty of the methodology is there- 
fore related to the theoretical estimate of the r-production 
ratio (Th/Eu) r . 

Unfortunately, the r-process remains the most compli- 
cated nucleosynthetic process to model from the astro- 
physics as well as nuclear physics point of view (for a 
review see Arnould & Takahashi 1999). On the nuclear 
physics side, the nuclear structure properties (such as the 
nuclear masses, deformation, . . . ) of thousands of nuclei 
located between the valley of /3-stability and the neutron 
drip line have to be known, as well as their interaction 
properties, i.e the (n, 7) and (~f,n) rates, a- and /3-decay 
half-lives and the fission probabilities. Despite much re- 
cent experimental effort, those quantities for most of the 
nuclei involved in the r-process remain unknown, so that 
they have to be extracted on theoretical grounds and are 
subject to the associated uncertainties. On top of these 



nuclear difficulties, the question of the astrophysical con- 
ditions under which the r-process can develop is far from 
being settled. The site(s) of the r-process is (are) not iden- 
tified yet, all the proposed scenarios facing serious prob- 
lems. For this reason, only parametric approaches, such 
as the so-called canonical model (Seeger et al. 1965) can 
be used to estimate the Th production. We use in the 
present study the multi-event model[] (Bouquelle et al. 
1996; Goriely & Arnould 1996) in which the best fit to the 
solar abundances is derived from a superposition of canon- 
ical events with the aid of an iterative inversion procedure. 
Compared with other treatments of the canonical model 
(e.g Pfeiffer et al. 1998), a major advantage of the multi- 
event approach is to provide an efficient tool for a system- 
atic study of the various uncertainties affecting the model 
(Goriely 1999). The iterative inversion method works in 
such a way that the modification of a given (nuclear or 
astrophysics) input in the r-process model leads to an au- 
tomatic renormalization of the thermodynamic conditions 
necessary to optimize the fit to the solar r-abundance dis- 
tribution. Therefore, the uncertainties affecting the input 
data of the parametric model, as well as their impact on 
the Th production can be studied systematically within 
the multi-event approach, as shown in the next section. 

Our standard calculation is performed under the fol- 
lowing thermodynamic conditions: 1.3 < T[10 9 K] < 1.7, 
10 22 < iV n [cm- 3 ] < 10 29 and 10 < n cap < 200 (where 
T is the temperature, N n the neutron density and n cap 
the number of neutrons captured per seed nucleus). Note 
that the r-process calculations are performed making use 
of the waiting point approximation, since under the ther- 
modynamic conditions considered here, an almost com- 
plete (n, 7) — (-f,n) equilibrium is established (Goriely & 
Arnould 1996). When not available experimentally, the 
nuclear data are taken from the ETFSI nuclear masses of 
Aboussir et al. (1995) and from the gross theory (GT2) of 
(3~ decay (and /3-delayed neutron emission) of Tachibana 
et al. (1990). In addition, a-decay and fission processes are 
also considered (before and after the neutron irradiation 
freeze-out). The fission processes include spontaneous, (3- 
delayed and neutron-induced fission, the probabilities of 
which are calculated according to the prescriptions of 
Kodoma & Takahashi (1975) with the ETFSI fission barri- 
ers (Mamdouh et al. 1998). The procedure used to fit the 
solar r-abundance distribution is similar to the one de- 
scribed in Bouquelle et al. (1996) though each isotope is 
now given a weight inversely proportional to the error af- 
fecting its solar r-abundance (Goriely, 1999). Since we are 
mainly concerned with Th cosmochronometry, no details 

1 Note that, in the case of r-processes responsible for ele- 
ments observed in ultra-metal-poor stars, the denomination 
"multi-event" does not refer to numerous astrophysical events, 
such as supernova explosions, but rather to numerous compo- 
nents of a given astrophysical event characterized by different 
thermodynamic conditions, for example in the different layers 
of a given supernova. 
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are given for the representative thermodynamic conditions 
required to fit the solar system r-abundance distribution 
(such details can be found in Goriely and Arnould, 1996). 

3. Uncertainties in the predicted Th abundance 

The r-process production of Th is obviously model de- 
pendent. However, for cosmochronological purposes, it is 
of fundamental importance to know to what extent the 
remaining uncertainties in the r-process modelling can af- 
fect the Th synthesis. From Eq. (|j), it can be seen that 
the Th abundance has to be determined within less than 
16% if we hope to predict the age of star within less than 
3 Gyr. A high accuracy has already been achieved obser- 
vationally with errors reduced to loge = 0.08 affecting the 
stellar age by about 3.7 Gyr. Unfortunately, other uncer- 
tainties still need to be solved. These mainly concern the 
r-process modelling, but before focussing on this subject, 
it is of interest to stress that the normalization to the Eu 
abundance is not free from uncertainties. As shown by 
Goriely (1999), even if the r-process models would be able 
to reproduce exactly the Eu solar r-abundance, this value 
is still uncertain by about 20% (the observed abundance 
in the ultra-metal-poor stars is known within 35%) lead- 
ing to an error in T» of about 6.6 Gyr. Note that with 
respect to normalization procedure, it might be safer to 
use the Ho abundance, since Ho is made of one stable iso- 
tope only and the s-contribution to its solar abundance is 
even smaller than for Eu, i.e the error bars on its solar 
abundance (about 13%) are smaller than in the Eu case. 
The additional uncertainties related to the predicted Eu r- 
abundance are neglected in the present study by normaliz- 
ing the calculated Th abundance to the solar r-abundance 
of Eu. The complete absence of correlation in the pro- 
duction of Th and Eu in the canonical approach of the 
r-process justifies this choice. 

The sensitivity of the calculated r-process abundances, 
and in particular of the Th abundance, to the different cru- 
cial inputs used in the multi-event model is now examined 
and the impact of such uncertainties on the age of the CS 
22892-052 star (Tf s ) is discussed. 

3.1. Sensitivity to astrophysics conditions 

Among the different thermodynamic parameters entering 
the canonical model, the most critical one affecting the 
Th synthesis is obviously the maximum number of neu- 
trons captured by the initial seed nuclei, n™p X , which de- 
fines the strong component of the r-process. In analogy 
with the s-process nucleosynthesis, we can define a main 
r-process component responsible for the production of all 
the elements up to the A = 195 peak and part of the Pb 
peak. This main component requires a value of n cap up 
to about 140. Till now, there is no constraint from real- 
istic models on the largest value that n cap can take, i.e 
in analogy with the s-process, on a strong r-componcnt 



responsible for the bulk production of Pb and Bi. Con- 
sidering canonical events with values of n cap > 140 would 
lead to the production of the Pb-peak elements, as well 
as Th, without affecting the synthesis of the lower-mass 
elements. To illustrate such a sensitivity, multi-event cal- 
culations are performed considering canonical events with 
a maximum number of neutrons captured of n™^ v x =140, 
145, 150 and 200 (Fig. |). An excellent fit is obtained for 
all isotopes with A < 204 and seen not to be affected at 
all by the change in the maximum value of con- 
sidered. On the contrary, increasing the above 140 
leads to an increase in the production of the Pb-peak, Th 
and U elements. The large uncertainties in the Pb-peak 
r- abundances cannot favour one or another fit. The com- 
plete absence of a strong r-component, i.e a maximum 
value of n cap = 140, leads to a negative age (when derived 
from Eq. Q) of the CS 22892-052 star, and can obviously be 
rejected. Including a strong r-component with n cap = 145, 
150 and 200 leads to a Th abundance that implies a star 
age T° s = 12.2, 22.9 and 28.9 Gyr, respectively. Since 
the fit is constrained by the upper value of the Pb and 
Bi abundances, increasing nj^f above 200 does not affect 
the upper value of the Th abundance. More precisely, no 
event with a value of n cap > 170 contributes to the fit to 
the solar system abundances with our adopted nuclear in- 
puts. In summary, any age below about 29 Gyr can thus be 
obtained just by adjusting the strength of such a strong 
r-process component, unless the s- or r-origin of the Pb 
and Bi can be determined with a greater accuracy. 

3.2. Sensitivity to nuclear physics input 

The most fundamental nuclear input to r-process models 
is well known to be the nuclear masses. Various mass mod- 
els are available, but for practical reasons we only consider 
here in addition to the ETFSI model, the ETFSI-Q model 
(Pearson et al. 1996) which takes into account the strong 
shell-quenching found in some microscopic calculations on 
highly neutron-rich nuclei, the popular FRDM model of 
Moller et al. (1995) and the recently-developed model of 
Duflo & Zuker (1995), hereafter DZ, based on a very differ- 
ent approach than the previously cited models and which 
has proven its remarkable ability to predict experimentally 
known masses. Many studies have compared the quality 
of these models and their differences in the prediction of 
masses far away from the valley of /3-stability. Their im- 
pact on the r-process nucleosynthesis has also been ana- 
lyzed in various papers (e.g Goriely & Arnould 1992), so 
that we will restrict ourselves to analyze their respective 
predictions of the Th abundance. 

Multi-event calculations are now performed making 
use of the 4 above-cited mass models. The resulting fits 
are shown in Fig. |^. The fits to the stable nuclei are of 
the same quality, and in particular in the Pb region, no 
major differences in the predicted r- abundances can be 
observed. In particular, it should be emphasized that no 
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Fig. 1. Comparison between the solar system r-abundances (Goriely 1999) and the distribution predicted by our 
standard multi-event superposition of events characterized by maximum values of =140 (dashed), 145 (dot-dash), 
150 (dotted) and 200 (full curve). The square corresponds to the Th abundance observed in CS 22892-052 (Sneden 
et al. 1996). The Bi and Th abundance are connected by a straight line to visualize the extrapolation predicted by 
the respective models. The vertical lines correspond to error bars in the solar system abundances taken from Goriely 
(1999). 



major deficiency in the fit is obtained in the pre-peak re- 
gions at A ~ 120 and A ~ 180 whatever mass model 
is used. Given the absence of realistic r-process models, 
there is obviously no reason to favour one or another 
mass formula on grounds of parametric fits to the solar r- 
abundance distribution, especially when dealing with the 
Th abundance predictions which exclusively depend on 
the r-process paths in the A > 232 region. The extrapo- 
lation to the Th abundance appears to be highly affected 
by the mass model used. The estimate of the star age 
amounts to T? s = 28.9, 23.8, 42.2 and 8.7 Gyr for the 
ETFSI, ETFSI-Q, FRDM and DZ models, respectively. 
Such differences are not surprising, since it is well known 
that the r-process paths for these mass models are sig- 
nificantly different, in particular in the Z > 82 region 
where the strength of the shell correction energy around 
the N = 184 shell closure can be very different. This is not 
the case for the ETFSI and ETFSI-Q models, because the 



shell quenching introduced in the ETFSI-Q model in the 
vicinity of the N — 184 shell closure is small. Therefore, 
the abundance predictions in the Pb and actinide regions, 
and consequently the stellar age predictions, are globally 
similar when making use of ETFSI or ETFSI-Q. Com- 
pared with the other models, the DZ formula is charac- 
terized by a steep slope of the mass parabola and a weak 
shell effect around N = 184, so that the progenitors re- 
sponsible for the final Pb abundance are found in a lower 
mass region by-passing partially Th and U. The Th abun- 
dance obtained with FRDM model is higher than with the 
ETFSI mass models, because of a more widely spread shell 
effect in the vicinity of the N < 184 shell closure affect- 
ing the r-process path down to N = 170. The abundance 
peak around N = 184 before freeze-out is consequently 
flattened to lower masses than in the ETFSI case and is 
less affected by fission processes after freeze-out. A higher 
Th abundance predicted with the FRDM model leads to a 
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Fig. 2. Same as Fig. || where the predictions are obtained with the ETFSI (full), ETFSI-Q (dashed), FRDM (dot-dash) 
and Duflo & Zuker (dotted) mass models. 



higher age estimate. The high sensitivity of the predicted 
Th abundance to the mass model will not be resolved be- 
fore improving our mass predictions in the heavy (Z > 82) 
neutron-rich region. 

Another important ingredient in the Th nucleosyn- 
thesis concerns the fission processes, i.e the spontaneous, 
neutron-induced and /3-delayed fission. Most of the r-pro- 
cess calculations do not include the fission processes at 
all or only partially. However, when dealing with Th cos- 
mochronometry, fission processes must be included in the 
most careful way in order to describe the competing pro- 
cesses responsible for the final Th abundance (namely a- 
decays, /3-decays and fissions) correctly. The recent large 
scale calculation of ETFSI fission barriers (Mamdouh et 
al. 1998) up to A = 295 is used to test the significance 
of fission on the Th cosmochronometry. The spontaneous, 
neutron-induced and /3-delayed fission probabilities are de- 
termined in the same way as in Kodoma & Takahashi 
(1975). It should, however, be stressed that when not avail- 
able experimentally, the spontaneous fission rates are de- 
rived from a new regression fit to experimental data based 
on our fission barrier predictions. Figure || shows the nu- 



clear regions where the different fission modes influence 
the r-process flows. Because of the strong ETFSI shell ef- 
fect on the fission barriers around N = 184, no fission 
recycling is found during the neutron irradiation, at least 
before crossing the N = 184 closure. Only r-process paths 
characterized by an astrophysical parameter S a ^ 2 MeV 
(for more details about the astrophysical parameter, see 
Goriely & Arnould 1992) are stopped by neutron-induced 
fission. Spontaneous fission can also affect such r-process 
paths before neutron freeze-out. /3-delayed fission is found 
to be of small importance compared with the other de- 
caying modes, even after freeze-out. On the contrary, the 
spontaneous fission is found to be faster than the /3-decay 
for almost all isobaric chains above A — 250 and crucial in 
estimating the final r-abundances in the Pb and actinide 
region. In the specific fits studied in the present paper, the 
fission fragments do not affect the low-mass abundance 
distribution. Obviously, all the above conclusions should 
be taken with care, because of the uncertainties remain- 
ing in the determination of the fission barriers and fission 
probabilities, which are to be studied in a forthcoming 
paper. 
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as well as the neutron drip line (dashed line). 



In order to quantify the impact of the fission processes 
on the Th cosmochronometry, a multi-event calculation 
is reiterated switching off all the fission processes. This 
numerical test just aims at illustrating the largest error 
possibly made when neglecting fission, but should not be 
regarded as a sensible test case for the Th prediction. It is 
found that the neglect of fission gives rise to an increase 
of the age of CS 22892-052 by 12.3 Gyr on grounds of 
the abundance distribution shown in Fig. ^. It can also be 
seen that when including fission processes, the fission frag- 
ments do not modify the global abundance distribution. 
Obviously, a complete and consistent treatment of the fis- 
sion processes (especially spontaneous fission) is required 
to build a reliable cosmochronometry on the actinides. 

Figure || also presents uncertainties associated with j3- 
decays and /3-delayed neutron emission by comparing the 
solar fits obtained with the GT2 model and the QRPA 
model of Moller et al. (1997). Both models have been ex- 
tensively used in previous works dedicated to the r-process 
nucleosynthesis, so that it is of interest to study their in- 
fluence on the Th cosmochronometry. If use is made of the 
QRPA model instead of the GT2 model, a reduction from 
28.9 Gyr down to 15.1 Gyr is obtained for the age of CS 
22892-052. Once again, it should be added that although 



the fit to the solar distribution obtained with the QRPA 
model is slightly worse than the one obtained with the 
GT2 approach, it cannot be rejected a priori, since other 
nuclear or astrophysics shortcomings of the model can be 
responsible for the observed discrepancies (for example in 
the A = 180 region). As stated previously, given our poor 
understanding of the r-process nucleosynthesis (especially 
of the astrophysical site) the quality of nuclear models 
should not be tested on astrophysics arguments like fits 
to the solar abundance distribution. 

3.3. Uncertainties in the solar r-abundance 

The uncertainties still affecting the s-process model are 
responsible for non-negligible imprecisions in the determi- 
nation of the residual r-abundances of the solar system 
content (Goriely 1999), in particular for the so-called s- 
dominant isotopes and Pb-Bi isotopes (see the large er- 
ror bars in Fig. [I]). The principal source of uncertainty in 
the solar r-abundance of Pb and Bi lies in our ignorance 
of the relative s- and r-contribution to their production. 
Both processes can produce Pb-Bi almost entirely, so that 
the r-contribution cannot be estimated in a reliable way on 
grounds of s-process calculations. This problematic aspect 
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Fig. 4. Same as Fig. [l] where the predictions are obtained using the GT2 /3-decay and /3-delayed neutron emission 
rates with (full) or without (dash-dot) fission processes. The dotted curve corresponds to the use of the QRPA /3-decay 
and /3-delayed neutron emission rates including fission. 



of the solar abundance splitting in the Pb region can only 
be resolved through realistic modelling of the s-process 
(or accurate abundance determination of Pb at the sur- 
face of ultra-metal-poor stars provided the assumption of 
the r-process universality be confirmed) . It should be kept 
in mind that the prediction of the Th and Pb-Bi abun- 
dances are strongly correlated, so that any uncertainty in 
the solar r-abundances of the Pb and Bi elements is trans- 
lated into the exponentially dependent uncertainty in the 
star age. Among the A > 200 r-isotopes, 204 Hg, 206 Pb 
and 209 Bi play an important role since their r-abundances 
are better determined than for their neighbours. The so- 
lar r-abundance of 204 Hg is indeed well determined, so 
that it seems logical to constrain the fit in such a way 
as to reproduce the 204 Hg solar r-abundance. As regards 
206 Pb and 209 Bi, they are characterized by a relatively 
well-determined abundance to which the Th production 
is directly correlated. However, reproducing the recom- 
mended solar abundance of 204 Hg, 206 Pb and 209 Bi si- 
multaneously appears to be impossible without strongly 
deteriorating the fit to the A = 195 peak. For this reason, 
we reiterate multi-event calculations in which the fitting 



procedure is constrained (with an extra statistical weight) 
to the solar abundance of A = 204 and 206 in one case 
and A = 204 and 209 in the other case, using two different 
mass models, namely the ETFSI and DZ models (Fig. ||). 
Although a negative age is obtained with the DZ masses 
when constraining the fit to 209 Bi, the other predicted 
ages are 30.4 Gyr for the ETFSI calculation constrained 
to 206 Pb and 10.5 Gyr in the two remaining cases. So, in 
addition to the large sensitivity of the stellar age to the 
mass models as studied in the previous section, the un- 
certainties in the solar r-abundances appear to affect the 
age determination by about 20 Gyr. As long as the s- or 
r-origin of the Pb and Bi solar abundance is not deter- 
mined with high accuracy, the Th cosmochronometry will 
not provide any reliable age estimate. 

4. The age of the stars 

So far, we have been evaluating the stellar ages on the 
basis of a fit to the solar r-abundance distribution. Such 
a procedure is obvioulsy based on our initial fundamen- 
tal assumption that the r-process site is unique in the 
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Fig. 5. Same as Fig. [j] when the fit is constrained on 204 Hg and 206 Pb solar abundances with the ETFSI (full) or 
DZ (dashed) mass models. The dash-dot and dotted lines correspond to the constrain on 204 Hg and 209 Bi abundances 
with ETFSI and DZ mass models, respectively. 



Galaxy. This basic assumption is a fundamental prerequi- 
site to build a Th cosmochronometry upon the abundance 
analysis of metal-poor stars at the present time, since, as 
explained above, a direct fit to the abundance distribution 
observed at the surface of metal-poor stars would present 
even larger uncertainties, because of the restricted num- 
ber of elements observed, the impossibility to distinguish 
isotopic ratios and the much smaller accuracy in the abun- 
dance data as compared with our solar system. However, 
in the specific case of Pb, the uncertainties still affecting 
the solar abundance might be greater than the ones found 
in the observation at the surface of HD 115444 or HD 
126238 (Sneden et al. 1998), so that the Pb abundance in 
such stars could probably be more constraining on the Th 
predictions than the solar value. Unfortunately, the Th 
abundance has not been determined yet in such stars, so 
that in this case, we limit the discussion on the relative 
stellar age. Compared with CS 22892-052, only a small 
number of r-elements are observed in HD 115444, but it is 
at the moment the only metal-poor star in which Pb and 
Th lines are detected simultaneously. All the physically 
sound (i.e leading to a positive age T^ s ) calculations pre- 



sented in the previous section are compared, in Fig. |[ 
with the elemental abundance distribution observed in 
HD 115444 and CS 22892-052. We only retain calcula- 
tions which provide a good fit to the observed Z > 55 
abundances, and in particular to the Pb abundance in the 
case of HD 115444. In order to achieve a good fit to all 
elements observed and to optimize the extrapolation to 
the Th region, a normalization of the abundance curves 
is done on the heaviest elements accurately observed, i.e 
Pt and Os for HD 115444 and CS 22892-052, respectively. 
In this case, the age of CS 22892-052 is found to lie in 
the 7 < T^fGyr] < 39 range, and a similar error range 
of about 32 Gyr is predicted for the age of HD 115444. 
The different predictions of the stellar age, as well as the 
elemental abundances of Eu, Pb, Bi, Th and U are sum- 
marized in Table |l|. The stellar age can be determined 
in two different ways. If normalized to the calculated Eu 
abundance, the age T 1 has the advantage of beeing free 
from uncertainties in the normalization procedures on the 
observed abundances, but is sensitive to the theoretical 
uncertainties made in the r-process predictions. In partic- 
ular, it is well known that the origin of the r-nuclides in 
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Fig. 6. Elemental abundance distributions observed in HD 115444 and CS 22892-052 (open squares) compared with 
multi-event r-abundance predictions (see text). 

where r(U) = 6.41 Gyr is the characteristic a-decay time- 
scale of U. As seen in Fig. |6| and Table 1, log(Th/U) r is 
found to lie within an 0.1 range, whatever theoretical in- 
puts are used in the r-process model. Such an accurate 
estimate is principally bound to the fact that Th and U 
are neighbour nuclei, and consequently their production 
ratio is not strongly affected by unreliable extrapolation 
procedures, but rather by local nuclear uncertainties, such 
as nuclear masses or fission processes in the actinide re- 
gion. 

From Eq. (0), it is found that a ±0.1 error on the ob- 
served or predicted production ratio of Th/U gives rise 
to a ±2.1 Gyr error on T*. A future simultaneous obser- 
vation of Th and U lines in ultra-metal-poor stars could 
therefore open the way to an accurate age determination 
in contrast to the still complicated Th cosmochronome- 
try based on Th lines only. Equation is shown in a 
graphical form in Fig. [?] for 3 possible values of the Th/U 
production ratio. Note, however, that a deeper analysis 
of the Th/U production ratio would be required before 



the A ~ 160 region is not easily explainable (e.g Meyer 
& Mullenax 1998). The impact of such errors remains to 
be estimated. On the contrary, normalizing the calculated 
Th abundance on the observed Eu avoids theoretical com- 
plications inherent in the origin of the A ~ 160 nuclides, 
but gives rise to additional errors of the order of ±0.1 dex 
on the abundances (as seen in Fig. ^), i.e of ±4.7 Gyr on 
the age T?. 

Although the present study suggests that at the mo- 
ment great care should be taken in estimating the age 
of the stars on the basis of observed Th abundance, new 
accurate observations of heavy r-elements could put the 
Th cosmochronometry on safer grounds, especially if Th 
and U lines could be observed accurately and simultane- 
ously in metal-poor stars, as already stressed by Arnould 
& Takahashi (1999). As a matter of fact, if the Th and U 
lines are available, an age estimate could be derived from 
the expression 
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Table 1. Ages T* (in Gyr) of CS 22892-052, elemental abundances (in loge) for Eu, Pb, Bi, Th, U and log(Th/U) 
predicted by the different theoretical calculations shown in Fig. || and analyzed in Sect. 3. T£ (T?) corresponds to the 
age calculated with a Th abundance normalized to the calculated (observed) Eu abundance. 



Comment (see Sect. 3) 


Eu 


Pb 


Bi 


Th 


U 


Th/U 


Ti 


rj-i'Z 


Standard (n™ 1 = 200, ETFSI, GT2) 


-0.79 


-0.20 


-0.76 


-1.00 


-1.16 


0.16 


22.9 


25.5 


n zi = 145 


-0.82 


-0.43 


-1.03 


-1.38 


-1.49 


0.11 


6.6 


8.2 


^cap 150 


-0.80 


-0.27 


-0.85 


-1.13 


-1.23 


0.10 


17.3 


19.5 


ETFSI-Q masses 


-0.79 


-0.26 


-0.86 


-1.10 


-1.32 


0.22 


18.2 


21.0 


FRDM masses 


-0.91 


-0.14 


-0.64 


-0.72 


-0.67 


-0.05 


41.6 


38.9 


DZ masses 


-0.75 


-0.19 


-0.85 


-1.40 


-1.44 


0.04 


2.4 


6.9 


ETFSI + constraint on A = 206 


-0.86 


-0.22 


-0.81 


-1.04 


-1.18 


0.14 


24.3 


23.8 


DZ + constraint on A — 206 


-0.75 


-0.20 


-0.85 


-1.36 


-1.38 


0.02 


4.3 


8.9 


ETFSI + constraint onA = 209 


-0.74 


-0.48 


-1.11 


-1.34 


-1.36 


0.02 


4.7 


9.6 
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T, [Gyr] 

Fig. 7. Relation between the observed Th/U ratio and 
the age of the star T* for 3 different estimates of the pro- 
duction ratio. 

rushing into an age determination. In particular, uncer- 
tainties in the fission processes have not been included in 
the present study, but could possibly affect the predicted 
Th/U ratio. 

5. Conclusions 

The present paper analyzes critically the impact of the 
remaining uncertainties in the existing r-process models 
on the determination of the age of the Galaxy derived on 
grounds of the Th cosmochrometry. Although the direct 
and accurate observation of the Th abundance in ultra- 
metal-poor stars provides a possible method to date the 
age of the star without calling for complex model of the 
chemical evolution of the Galaxy, the Th cosmochronom- 
etry remains affected by all the difficulties associated with 
our poor understanding of the r-process nucleosynthesis. 
Even if we disregard the difficulty associated with the 
fundamental assumption relative to the unicity of the r- 
process site in the Galaxy, it is impossible at the present 
stage reliably to determine the age of a ultra-metal-poor 
star from parametric models of the r-process. The major 
difficulties lie in the very s- or r-origin of the Pb and Bi 



isotopes, in addition to the unavoidable problem related 
to the prediction of nuclear data (principally masses) and 
the still unknown thermodynamic conditions in which the 
r-process takes place (in particular concerning the maxi- 
mum neutron irradiations that can be reached). New accu- 
rate observations including Pb and Th lines in ultra-metal 
poor stars could shed light on the ability of the r-process 
to produce the Pb-peak elements, and consequently con- 
strain the Th cosmochronometry. In particular, the Th 
cosmochronometry could be put on safer grounds, espe- 
cially if Th and U lines could be observed accurately and 
simultaneously in ultra-metal-poor stars. As regards nu- 
clear uncertainties, only theoretical as well as experimen- 
tal effort to come can help us achieve the required accuracy 
to pretend to derive the age of the Galaxy on the basis of 
the actinides cosmochronometry. 

Acknowledgements. The authors are grateful to M. Arnould 
for helpful discussions 

References 

Aboussir Y., Pearson J.M., Dutta A.K., Tondeur F., 1995, At. 

Data Nucl. Data Tables 61, 127 
Arnould M., Takahashi K., 1999, Rep. Prog. Phys. 62, 393 
Bouquelle V., Cerf N., Arnould M., Tachibana T., Goriely S., 

1996, A&A 305, 1005 
Cowan J.J., Mc William A., Sneden C, Burris D.L., 1997, ApJ 

480, 246 

Duflo J., Zuker A. P., 1995, Phys. Rev C52, 23 
Goriely S., Arnould M., 1996, A&A 312, 327 
Goriely S., Arnould M., 1997, A&A 322, L29 
Goriely S., 1999, A&A 342, 881 

Kodoma T., Takahashi K., 1975, Nucl. Phys. A239, 489 
Mamdouh A., Pearson J.M., Rayet M., Tondeur F., 1998, Nucl. 

Phys. A644, 389 
Meyer B.S., Mullenax D.J. 1998, Nuclei in the Cosmos dol V, 

eds N. Prantzos, S. Harrisopoulos (Gif-sur-Yvette: Editions 

Frontieres) p. 289 
Moller P., Nix J.R., Myers W.D., Swiatecki W.J., 1995, At. 

Data Nucl. Data Tables 59, 131 
Moller P., Nix J.R., Kratz K.-L., 1997, At. Data Nucl. Data 

Tables 66, 185 



S. Goriely & B. Clerbaux: Uncertainties in the Th cosmochronometry 



Pearson J.M., Nayak R.C., Goriely S., 1996, Phys. Lett. B 387, 
455 

Pfeiffer B., Kratz K.-L., Thielemann F.-K. et al., 1998, In Proc. 

of the 9th Workshop on Nuclear Astrophysics, eds W. Hille- 

brandt, E. Miiller (MPA, Munich) p. 168 
Seeger P.A., Fowler W.A., Clayton D.D., 1965, ApJS 11, 121 
Sneden C, Mc William A., Preston G.W., et al., 1996, ApJ 

467, 819 

Sneden C, Cowan J. J., Burris D.L., Truran J.W., 1998, ApJ 
496, 235 

Tachibana T., Yamada M., Yoshida Y., 1990, Progr. Theor. 
Phys. 84, 641 



